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Abstract 



We compute the hadronic light-by-light scattering contributions to the 
muon anomalous magnetic moment, a" (had), in chiral perturbation theory 
that are enhanced by large logarithms and a factor of Nq- They depend on a 
low-energy constant entering pseudoscalar meson decay into a charged lepton 
pair. The uncertainty introduced by this constant is ±60 x 1CP 11 , which 
is comparable in magnitude to the present uncertainty entering the leading- 
order vacuum polarization contributions to the anomalous moment. It may be 
reduced to some extent through an improved measurement of the ir° — > e + e~ 
branching ratio. However, the dependence of a" (had) on non-logarithmically 
enhanced effects cannot be constrained except through the measurement of 
the anomalous moment itself. The extraction of information on new physics 
would require a future experimental value for the anomalous moment differing 
significantly from the 2001 result reported by the E821 collaboration. 
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The recently reported measurement of the muon anomalous magnetic moment a M by the 
E821 collaboration |l| has generated considerable excitement about possible evidence for new 
physics. The interpretation of the result, however, depends in part on a reliable treatment 
of hadronic contributions to a M which arise at two- and three-loop order. While the vacuum 
polarization contribution can be constrained by e + e~ experiments and r decays and appears 
to be under adequate theoretical control ||, recent analysis [|3],f§] of the hadronic light-by-light 
scattering contribution a^ i (had) have uncovered a sign error in previous calculations of 
the dominant, pseudoscalar pole term. The resulting sign change reduces the 2.6a deviation 
of a M from the Standard Model prediction reported in JTJ by one standard deviation, thereby 



modifying considerably the original interpretation of the result [|K||ri[ . 

The commonly quoted values for a" (had) (after incorporating the corrected overall sign) 
rely on model treatments of the off shell 7r°7*7*, 777*7* and 7/7*7* interactions. While the 
amplitude for pseudscalar decay into two real photons is dictated by the chiral anomaly, the 
off-shell amplitudes relevant for a" (had) are affected by non-perturbative strong interactions 
whose effects cannot yet be computed with sufficient precision from first principles in QCD. 
Similarly, the contributions from other hadronic intermediate states besides the ir°, rj, and 77' 
cannot be computed reliably at present. The analysis of a" (had) falls naturally under the 
purview of chiral perturbation theory (%PT), which provides systematic, model-independent 
framework for parameterizing presently incalculable hadronic effects. The static quantity 
a" (had) has an expansion in powers of p/A, where p is a small mass of order or m n and 
A is a hadronic scale, typically taken to be ~ AttF^ ~ lGeVQ. The coefficients appearing in 
the expansion depend in part on a priori unknown "low-energy constants" (LEC's), which 
parameterize the effects of non-perturbative short distance physics. In principle, the LEC's 
may be determined from an appropriate set of experimental measurements. 

In this Letter, we perform a x?T calculation of a" (had) including all the logarithmically 
enhanced contributions to a" (had) that arise at order Nca 3 p 2 /A 2 , where Nq is the number 
of quark colors. We identify the dependence of a" (had) on the large logarithms of A/p 
as well as on the relevant LEC's. The result for the large In 2 term - which is determined 
entirely by gauge-invariance and the chiral anomaly - was first given in Ref. and was 
used to uncover the sign error in previous model calculations. However, only part of the 
large In term is fixed by symmetry considerations. It receives an additional contribution 
involving x, a LEC entering the rate for pseudoscalar decay into leptons. The x-dependent 
piece was also computed in Ref. [|J. An analytic calculation of the x~ independent large In 
term was performed by the authors of Ref. ||12|| , who used a model for the off-shell P r y* r j* 



form factor to regulate the two- loop amplitude and assumed m M was almost equal to m^. 
The results of the latter calculation also contain a model prediction for the non-logarithmic 
0(N c a 3 p 2 /A 2 ) contribution to a" (had). 

In what follows, we carry out a consistent xPT treatment of a" (had), providing the first 
complete, model-independent analysis of the terms enhanced by large logarithms through 
0(Nc<y 3 p 2 /A 2 ). The sum of these terms is known, since existing measurements for 77 — > 
fi + H~ and 7T° — ► e + e~ branching ratios fix the value of \. However, the uncertainty in 
a" (had) from the error in x is significant: ±60 x 10 -11 , which is roughly the same size as 



1 We treat niu and to be of the same order and take both to be small compared with A. 
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the present theoretical uncertainty in the leading-order vacuum polarization contributions to 
a^. In principle, an improved measurement of the 7r° — > e + e~ branching ratio could reduce 
this source of uncertainty. 

A more serious consideration involves the contributions to a" (had) that are beyond 
the order at which we compute. These include effects of order 0(a 3 p 2 /A 2 ) which are not 
enhanced by a factor of Nc and order 0(Ncce 3 p 2 /A 2 ) terms that are not enhanced by large 
logarithms. We parameterize all these effects in terms of C . For simplicity we will refer to C 
as a low energy constant even though it has nonanalytic dependence on the quark masses. 
At present, C cannot be determined in a model-independent way without reliance on the 
measurement of a M itself. Even the sign of the C*-dependent term cannot presently be fixed 
in an model-independent manner. The presence of this LEC renders the interpretation of 
in terms of new physics problematic, since the size of the (7-dependent contribution could be 
as large as the present experimental error in a M Below we discuss the conditions under 
which one might still be able to extract information on new physics from an a M measurement. 

To leading order in Nc, the lowest-order chiral contributions to a" (had), enhanced by 
large logarithms, arise from the two- and one-loop graphs of Fig 1. Taking and m M 
as being of 0(p), the leading, large logarithmic contributions arise at order Nca 3 p 2 /A 2 . 
The two-loop graphs (Fig. la) contain an overall, superficial cubic divergence as well as 
a linearly-divergent one-loop subgraph involving two photons and a muon line. The latter 
must be regulated by adding the appropriate one-loop counterterm (ct) (Fig. lb). The 
one-loop graphs also contain an insertion of x(/ i )- The sum of these graphs contains a 
residual divergence, which must be removed by the appropriate magnetic moment ct (Fig. 
lc). Associated with this ct is a finite piece which, as discussed above, can only be fixed in a 
model-independent way by the measurement of a M itself. Additional contributions also arise 
from the graphs such as those appearing in Fig. 2. Although subdominant in Nc counting 
the contribution of the charged pion loop in 2b (and related diagrams) is leading order in 
p/A. The order 0(a 3 N c ) term arising from the three loop graphs with a charged pion loop 
- which we denote by a"(had)i. . - has been computed in Ref. [|15| ,|5ll- The result is finite 



and contains no large logarithms. We will include this contribution when we compare the 
theoretical prediction for a" (had) with experiment. 

At order a 3 p 2 /A 2 there will be contributions to a" (had) from higher dimension operators 
inserted at the vertices in the charged pion loop graph of Fig. 2 and from four-loop graphs 
containing an additional hadronic loop. Some of these should contain large logarithms. 
However, these are suppressed by a factor of N c compared to the logarithmically enhanced 
pieces that we compute. As remarked earlier we absorb these and many other subdominant 
contributions (e.g., Fig. 2a) into C and do not discuss them further here. 

As inputs for amplitudes of Fig. 1, we require the Wess-Zumino-Witten -P77 interaction 
Lagrangian 



as well as the leading-order operator contributing to the decays P — > £ + £ [[14| : 

*W = ! ^^/i7 A 75/i[xiTr(g 2 St J D A S - Q 2 D x ^Il) (2) 
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where 



+X 2 Tr (Q^QDj: - QD x tfQZ 



(3) 



(4) 



Here, i 7 ^ = 92.4 MeV is the pion decay constant, A a denote the Gell-Mann SU(3) matrices, 
and D\ is the covariant derivative. Note that in contrast to the conventions of Ref. we 
have made the iVc-dependence of the LEC's \i explicit for the sake of clarity. 

In computing the loop amplitudes involving these operators, it is important to employ 
a regulator which maintains the consistent power-counting of the chiral expansion. To that 
end, we employ dimensional regularization, where we continue only momenta (and not Dirac 
matrices) into d = 4 — 2e dimensions. The relation between bare and renormalized couplings 
is, 



xl + xl 
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Using Eq. fl5|) and adding the amplitudes for Fig. la,b, we obtain the divergent part of the 
two-loop amplitude 
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where q a and e@ are the photon momentum and polarization, respectively. We remove this 
divergence using a magnetic moment counterterm. The bare coupling Co and renormalized 
coupling C(//) are related by 
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The light-by-light contribution to the anomalous magnetic moment a" (had) is a physical 
quantity and has no dependence on the subtraction point \i. The //-dependence of the 
diagrams cancels that of the couplings C{p) and x(/-0- To obtain the /i-dependence of the 
couplings we require that the bare Green's functions corresponding to the sum of Figs, la-c 
and the to the Pi + i~ one-loop subgraphs, respectively, be independent of the subtraction 
scale. Doing so leads to a coupled set of renormalization group equations for an d 



dx 
d\x 
dC 
^ dfi 
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(9) 
(10) 



The solution is 
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X(n) = 121n(/i //i) + x(A*o) (11) 
C(/i) = 61n 2 (/i / / i) + [x(/io) + 3]ln( / io//i) + C( /Uo ) • (12) 

At a scale /io = A ~ 1 GeV, the constants C(fio) and xif^o) contain no large logarithms 
of the form ln fc (A/p) (k = 1,2) where p is around m M or m w . For /i of OQo), however, the 
Feynman diagrams contain no such large logarithms, and they live entirely in C(/x) and 
x(aO- Hence, the resulting expression for a" (had) is, in the MS scheme, 

a" (had) = a"(had)i. Q . 




7(0 + \ + ix(A) 



(13) 



where // is of order p and could be set equal to either m M or m^. Recall that A ~ 4:irF n ~ 1 
GeV. The function f(r), with r = m^/m 2 , arises from the one loop diagram with a coupling 
proportional to x(/-0 (Fig- lb) and is given byQ 
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/(r)=ln(j^J +-r 2 lnr--(2r + 13) + -(2 + r)^/r(4-r)cos^ J . (14) 

Note that we have absorbed all the remaining terms, including those proportional to C(A) 
and x(A) not enhanced by large logarithms, into C . 

The logarithmically enhanced hadronic light-by- light contributions to a M are renormal- 
ization scheme-independent. However, the values of x(A), f(r), and the constant appearing 
in the renormalization group equation for C (leading to the factor of 1/2 in the In term in 
Eq. (|I^) ) depend on one's choice of scheme^. This scheme- dependence cancels in the sum 
of their contributions to a" (had). For our calculations, we adopted a scheme in which the 
loop integrals were evaluated in <i-dimensions with d > 4 (corresponding to e < 0), while 
the Dirac matrices and photon polarization indices were taken as four-dimensional. For 
this choice we have tj^Tt^^^) = 16 instead of Ad. Moreover, the value of x(A) in this 



scheme is the same as that in |T4| but is four less than the xiA) used in | I7|| . An alternative, 
equally convenient scheme is again to treat the Dirac matrices, epsilon tensors, and photon 
polarizations as four dimensional; take d < 4 (e > 0); and rewrite (|3|) as 



K V>l»ll3lnV Vadx U° + -^VJ +■■■ > (15) 
where the metric tensor i] a! 3 is <i-dimensional. In this scheme the value of x(A) is still the 



. N c a 2 (xi + 
48tt 2 F^ 6 



same as m 



|14B , but f(r) — > f(r) + 3/2 and the —3 in the renormalization group equation for 
C becomes —12. Notice that the total value for the logarithmically enhanced contribution 
to a" (had) is unchanged. 



2 Since we only compute the terms enhanced by large logarithms to get / we replace x(a*) by 
121n(A//x). 

3 We will discuss the issue of scheme dependence and give more details of the calculation in a 
future publication. 
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As a check on the result in Eq. (jig), one may compute the one- and two- loop amplitudes 
with the insertion of x(A) in Fig. lb and C(A) in Fig. lc. In this case, all of the large 
logarithms arise from the Feynman amplitudes and not from the operator coefficients. Using 
an explicit calculation, we have verified in the limit m n — > that this procedure exactly 
reproduces the expression in Eq. ([13]). We note that the In term and the term proportional 
to x agree with the expression in Ref. ||. 

Chiral perturbation theory can be used for the r\ — > amplitude and the LEC 



x(A) can be deduced from the measured i] — > yU + /i~ branching ratio ||16|| . This yields (17 
x(lGeV) = —14^5 or —39^4, where we have scaled the results in Ref. [[H]] from A = m p up 
to A = 1 GeV and subtracted four. Note that because the rj — ► fi + n~ branching ratio is a 
quadratic function of %, two different values for this LEC may be extracted from experiment. 

Our calculation of the logarithmically enhanced contributions to a" (had) only involved 
the use of chiral SU(2)l x SU(2)r while this extraction of the LEC involves the use of 
chiral SU(3)l x SU(3)r. Since one expects chiral perturbation theory to work better in the 
case where only the pions are treated as light it is desirable to have an extraction of x that 
only relies on chiral SU(2)i x SU(2)r. This can be done using the measured 7r° — ► e + e~ 
branching ratio which yields [|17|] ^(lGeV) = — 29-11 or ^^-It- Unfortunately, the errors 
on the extracted x(lGeV) are very large in this case. A more precise determination of the 
7T° — > e + e~ branching ratio could reduce the theoretical uncertainty in x(lGeV). 

Model calculations for a" (had) differ from our analysis typically through insertion of 
form factors at the P7*7* vertices obtained from the WZW interaction in Eq. (0). For ex- 
ample, one widely-followed model employs form factors based on a vector meson dominance 
picture. This approach - known as resonance saturation - may also be used to obtain \i 
giving JHJ x(lGeV) res sa ^ ~ — 17 . In a similar vein, one may analyze this LEC at leading 
order in Nc, where it depends on a sum over an infinite tower of vector meson resonances 
|18| . Using a model-dependent form factor for the sum over vector resonances that at short 
distances is consistent with the properties of QCD gives x(lGeV) iVo: res sa ^ = — 16 ± 5. 
This result provides some quantitative support for phenomenological models since it is close 
to the value x(lGeV) ~ —14 obtained from experiment. However, one would not want to 
draw conclusions about the validity of the Standard Model using such a model-dependent 
approach. 

Using x(lGeV) = —14^5 as input, setting \x = m M , and adding the large In 2 and In terms 
in Eq. (|T3"D, we obtain 

a"(had) log =(57^)xl0- n . (16) 

We emphasize that inclusion of both the —f(r) + 1/2 and x/6 parts of the In (A///) term 
is crucial to obtaining an accurate numerical result for a"(had)i og . If, for example, one 
were to keep only the dependence on x(A), one would find substantial cancelations between 
the In 2 and In contributions. The presence of the calculable —f(r) + 1/2 term, however, 
substantially mitigates these cancelations. 

We observe that the central value in Eq. ([II]) is roughly a factor of two larger than 
obtained in model calculations for the 7r° contribution, and that the uncertainty is about 
a third the size of the present experimental error in a M JTJ] . After the full E821 data set is 
analyzed, the uncertainty in Eq. ( |I6|) will be comparable to the anticipated experimental 
error. As noted above, improved measurements of the tt° — > e + e~ branching ratio could 
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reduce the theoretical uncertainty in the large logarithmic contributions to a" (had). Using 
the other value of x obtained from the 77 — > branching ratio, x(lGeV) = —39^4, leads 

to a"(had)i og = (— 190^5q) x 1CT 11 . Although there exists a strong theoretical prejudice 
in favor of the first solution [Eq. ([16])] based on both the resonance saturation model for x 
and a" (had) as well as consistency between the values of x obtained from the r\ — > ji + \f~ 
and 7T — > e + e~ branching ratios, one cannot rule out the second value for a"(had)i og . 
Adding in the 0(N^a 3 ) charged pion loop contribution j5,|15|, a"(had)i. . = —44.6 x 



10 11 , to that in Eq. ( |T5] ) gives the following x?T expression for a" (had), 

'(had) = (l3±|g +31(7) x lO" 11 . (17) 



The largest uncertainty in the expression for a" (had) above arises from the subdomi- 
nant terms that have not been calculated and are parameterized by the LEC C. As noted 
above, this constant includes the effects of non-logarithmically enhanced two-loop contri- 
butions, heavy mesons such as the 77 and 77' which have been integrated out, and other 
non-perturbative dynamics. On general grounds, one could expect its natural size to be of 
order unity. A comparison of Eq.(|17]) with the results of model calculations is consistent 
with this expectation. For example, the model calculation of Ref. |TI[ corresponds roughly 



to C ~ 1. Rigorously speaking, however, the precise value - as well as the sign - of C 
is unknown. An uncertainty AC = ±1 corresponds to Aa"(had) = ±31 x 10 , which 
is roughly one fifth of a standard deviation for the published Brookhaven measurement [[!]. 
One should not, however, treat this as an estimate of the theoretical uncertainty in a" (had). 
A value of C equal to +3 or —3, for example, would not be unusual. 

Alternatively, one may use the experimental result for a M to determine C. To that end, 
we use the up-dated results for hadronic vacuum polarization contributions [H, the QED 



and electroweak loop contributions in Ref. (T^j and the value for a"(had) c hi ra i given in Eq. 
(|I~7j). ^From the E821 result for we obtain 

(7 = 7 + 5 + 3+2 , (18) 

where the first uncertainty arises from the experimental error in a M , the second corresponds 
to the theoretical QED, electroweak, and hadronic vacuum polarization errors, and the final 
uncertainty arises from the error in x- I n the future, the first uncertainty will be considerably 
reduced upon complete analysis of the full E821 data set. The value of C is consistent with 
unity, though it could be considerably larger, given the other experimental and theoretical 
inputs into Eq. (|l~8|) . Using the second solution for x and a" (had) gives C = 16 + 5 + 3 ±2. 

At present there is no indication that the hadronic LEC (7 differs substantially from its 
natural size and, thus, no reason to discern effects of new physics, such as loops containing 
super symmetric particles ]nj, from the a M result. In principle, a systematic calculation of 
some of the effects arising at higher order - such as terms of 0(a 3 N(jp 2 /A 2 ) enhanced by 
large logarithms - could modify this conclusion PU| . Similarly, should the full E821 data 



imply a value for (7 which differs significantly from ±1 (e.g., by an order of magnitude), 
one might argue that there is evidence of new physics. Such a conclusion would presumably 
require considerable disagreement between the published E821 result [|l| and the analysis of 
the full data set. The most convincing analysis, however, would rely on a first principles 
QCD calculation of a" (had), a prospect which seems to lie well into the future. 
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FIGURES 




(a) 




- X . — 

(o) 

FIG. 1. Hadronic light-by- light contributions to muon anomalous magnetic moment, a" (had), 
obtained from insertion of the WZW interaction at the 7r°77 vertices. The x in (b) indicates the 
insertion of the the low-energy constant x P ms its counterterm. The x in (c) denotes the magnetic 
moment coupling C plus its counter term. The solid, dashed, and wavey lines denote the fi, 7r°, 
and 7, respectively. 




h\ + 

(b) 

FIG. 2. Some hadronic light-by- light contributions to muon anomalous magnetic moment that 
are not enhanced by large logarithms. 
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